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INTRODUCTION 


Cased  explosive  munitions  can  respond  to  ballistic  impact  in  a 
variety  of  ways:  deformation,  ignition,  burning,  transition  from  burn¬ 
ing  to  a  violent  reaction  (after  a  variable  time  delay),  detonation,  or 
possibly,  no  reaction  at  all.  Because  of  the  variety  of  impact  condi¬ 
tions  and  system  responses,  it  has  been  difficult  to  quantify  the 
important  parameters  leading  to  reaction  and  the  prediction  of  the 
ultimate  system  response. 

Thermal  and  venting  situations  play  key  roles  in  determining  the 
fate  of  an  impacted  munition.  Sewell  and  Kinney  (Reference  1)  have 
studied  in  detail  the  phenomenon  of  venting.  They  later  developed  an 
equation  (Reference  2)  to  predict  the  critical  vent  size  necessary  to 
keep  a  burning  system  from  pressurizing.  This  equation  was  applied  to 
a  generic  explosive,  with  representative  constants  in  the  equation. 
The  results  of  this  calculation  provided  the  foundation  for  a  set  of 
confirming  experiments  for  a  real  explosive  (Composition  B).  The 
latest  series  of  tests  (Series  III)  is  reported  in  the  following 
sections. 

The  purpose  of  this  work  is  to  address  the  following  questions. 

1.  What  is  the  relationship  between  the  critical  vent  area  and 
the  burning  surface  area  of  the  explosive? 

2.  What  role  does  the  temperature  of  the  explosive  play  in  the 
calculation  of  the  critical  vent  area? 

3.  What  properties  of  the  vented,  burning  explosive  system  lead 
to  delayed  violent  reactions? 

4.  What  constants  in  the  generic  venting  equation  (Reference  2) 
need  to  be  modified  for  the  explosive,  Composition  B? 


VENTED  EXPLOSION  CHAMBER  (VEC) 
TEST  APPARATUS  AND  MATERIALS 


Explosive  samples  were  contained  in  1014  seamless  steel  cylinders 
that  were  5-1/2  inches  tall  with  a  1/2-inch-thick  wall  and  a  4-inch 
inside  diameter.  The  bottom  plate  was  electron-beam  welded  to  the 
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cylinder  to  provide  a  continuous,  integral  "cup-shaped"  container; 
while  the  top  plate,  containing  the  vent  hole,  was  bolted  at 
eight  places  to  the  cylinder  with  1-inch  long,  1/4  x  20  cap  screws. 
Static  pressure  to  separate  the  lid  from  the  container  was  about 
1400  psi.  (See  Figure  1.) 


i  a  *  «•  J 


•1  *  —  < '  '  V  i  <  ■ 


y«;. 


FIGURE  1.  VEC  Test  Apparatus. 


In  the  most  recent  Series  III  tests  (reported  here),  the  top  of 
the  apparatus  was  modified  to  include  an  ignition  feed-through  and  a 
Bell  and  Howell  Model  4-313A,  0-100  psi,  water-cooled  pressure  trans¬ 
ducer.  (See  Figure  2.) 

The  empty  cylinders  were  sent  to  the  Naval  Weapons  Station, 
Yorktown,  VA,  for  lining  and  explosive  loading.  The  1/8-inch-thick 
liner  was  a  polyethylene/polypropylene  material  used  in  the  NAVSEA 
Explosives  Advanced  Development  Program.  The  explosive  load  was 
nominally  2  pounds  of  cast  Composition  B  (density  =  1.70  g/cm3).  The 
explosive  filled  the  cylinder  to  a  height  of  3  inches  above  the 
1/8-inch-thick  liner,  leaving  a  1-7/8-inch  air  space  between  it  and 
the  lid  containing  the  vent  hole. 


FIGURE  2.  Series  III  Test  Apparatus. 


The  various  vent  holes  were  precisely  drilled  through  the  center 
of  the  1/2-inch-thick  x  5-inch-diameter  circular  lids.  The  vent  hoies 
functioned  as  square-edged  orifices  subject  to  compressible  flow. 

Reproducible  ignition  of  Composition  B  explosive  at  atmospheric 
pressure  proved  to  be  a  very  difficult  task.  Ignition  had  to  be 
accomplished  without  initially  overpressurizing  the  vented  container, 
which  would  lead  to  destruction  of  the  test.  Several  "gasless"  methods 
were  tried.  An  electrically  heated  Constantan  ignition  wire  ribbon, 
3/lb-inch  wide  by  =4-1/2  inches  long,  was  tried  without  success.  Next, 
a  combination  of  thermite  powder  and  thermite  ignition  mixture  was 
tried.  This  gave  a  short-duration,  high-heat  pulse  that  sometimes 
ignited  the  Composition  B,  but  more  often  behaved  like  a  dust  explo¬ 
sion,  which  rapidly  pressurized  the  system.  A  reproducible  ignition 
system  was  finally  developed  which  consisted  of  0 . 020-inch-thick  X 
3/4-inch-wide  Pyrofuse  (a  sandwich  of  palladium  and  aluminum,  which, 
when  heated  to  650°C,  alloys  to  PdA£2  at  2800°C).  This  was  connected 
(from  the  feedthrough  on  one  side  of  the  lid  to  the  metal  case  on  the 
other)  to  a  variable  arc  welding  power  supply.  Since  Pyrofuse  has  low 
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resistance,  a  high  current  of  200  A  was  necessary  to  heat  it  to  reac¬ 
tion  temperature.  Ignition  of  the  Pyrofuse  and  Composition  B  typically 
occurred  between  5  and  8  seconds  after  application  of  the  power  to  the 
Pyrofuse.  This  ignited  a  rectangular  area  on  the  surface  of  the  explo¬ 
sive  that  rapidly  spread  over  the  exposed  surface  of  the  Composition  B.  < 

The  initial  temperature  of  the  explosive  was  recorded.  For  low- 
and  high-temperature  runs,  rounds  were  conditioned  for  24  hours  at  the 
desired  temperature,  then  placed  in  a  well-insulated  container  for 
test  setup  and  actual  run. 


Pressure  measurements  versus  time  were  recorded  on  magnetic  tape 
for  later  oscillographic  reproduction.  Time  zero  was  set  at  the  igni¬ 
tion  of  the  Pyrofuse.  Real-time  video  tape  was  also  employed  to  aid  in 
determining  onset  of  reaction,  ultimate  reaction,  and  completion  of  the 
experiment.  Pre-  and  post-still  photography  were  also  used.  (Test 
parameters  are  listed  in  Table  1.) 


TABLE  1 .  Test  Parameters . 


Test  no. 

Serial  no. 

Vent 

diameter, 

in. 

Vent 
area , 
in2 

Explosive 

temperature, 

°C 

VEC  1 

8 

0.234 

0.0430 

15 

VEC  2 

9 

0.191 

0.0287 

15 

VEC  3 

6 

0.221 

0.0384 

0 

VEC  4 

5 

0.266 

0.0556 

61 

VEC  5 

7 

0.219 

0.0376 

63 

RESULTS 


The  experimental  data  are  presented  in  Appendix  A.  A  short  sum¬ 
mary  of  results  is  presented  in  Table  2. 
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TABLE  2.  Results  of  Series  III  VEC  Tests. 


Test  no. 

Burn  time,a 

s 

Reaction 

VEC  1 

269.0 

Slow  burn,  with  several  small  pressure  oscilla¬ 
tions  (P  >  300  psi  at  269  seconds) . 

(v  =  0.24  mm/s . )b 

VEC  2 

7.6 

Slow  burn  until  first  pressure  oscillation  at 

5.6  seconds.  (P  >  330  psi  at  7.6  seconds.) 

max 

VEC  3 

4.08 

Slow  burn  for  2  seconds,  then  slow  runup  until 
3.9  seconds.  Four  oscillations  in  pressure, 
then  pressure  rupture  (debris  scattered 

43  feet).  (Anomalous--later  examination  of 
high-speed  film  indicated  bottom  rupture — 
explosive  separated  from  liner.) 

VEC  4 

384.0 

Quiescent  burn  for  252  seconds.  Large  pressure 
buildup  at  267  seconds  followed  by  "hissing" 
and  13-inch  flame  out  the  vent  hole  from  276  to 
291  seconds.  Black  smoke  at  303  seconds. 

Reaction  completed  at  384  seconds. 

(v  -  0.17  mm/s . )^ 

VEC  5 

4.4 

Steady  pressure-rise  with  no  oscillations 
apparent.  At  4.4  seconds,  top  blew  off 
scattering  much  unreacted,  fractured  explosive 
over  a  wide  area. 

aAfter  ignition. 


v  is  the  average  linear  burning  velocity  of  the  explosive. 


DISCUSSION 


THEORY  OF  VENTED  BURNING 
Pressure  Rise  Rate 


Kinney  and  Sewell  (Reference  1)  determined,  from  interior  bal¬ 
listics,  the  rate  of  pressure  rise  from  combustion  of  an  energetic 
material.  The  basic  form  is  given  in  Equation  (1)  below: 


dp  _  _ B  dn 

dt  V  dt 


(1) 
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where  dn/dt  is  the  time  rate  of  change  of  the  number  of  moles  of 
product  gases.  This  equation  may  be  replaced  with  one  in  which  the 
variables  are  more  easily  measurable.  Thus, 


_ 

dt 


fi _ 9! _ 

M  (A  -  BTq) 


V 


(2) 


where 

-5  3 

R  =  molar  gas  constant  =  8.314  X  10  bar  -  m  /mol  -  K 

3 

V  =  volume,  m 

T„  =  flame  temperature,  K 
o 

M  =  formula  mass  product  gas,  kg/mol 

3 

p  =  density  of  explosive,  kg/m 
Tq  =  bulk  temperature  of  explosive,  K 
a,A,B  =  energetic  material  constants  (see  below) 

3 

Sg  =  burn  surface  area,  m 
P  =  absolute  pressure,  bars 


The  term  [a/ (A  -  BT  )]  represents  the  variation  in  burning  rate  with 
bulk  explosive  temperature.  The  values  of  the  generic  constants  were 
chosen  after  examining  the  work  of  Johansson  and  Persson  (Reference  3) , 
in  which  they  found  that  the  reciprocal  burning  rate  as  a  function 
of  temperature  for  porous  RDX  (p  =0.9  g/m3)  at  ambient  temperature 
was  about  0.5  mm/s.  However,  extrapolated  strand  burner  linear  burn¬ 
ing  rates  (Reference  4)  showed  a  lower  rate  of  burning  for  cast 
Composition  B  of  0.13-0.15  mm/s  at  ambient  temperature  and  pressure. 
The  infinite  burn  rate  intercept  was  estimated  as  510  K  (237°C) 
(Reference  3).  From  differential  thermal  gravimetry,  differential 
thermal  analysis,  and  thermal  gravimetric  analysis  experiments  at 
the  Naval  Weapons  Center  (NVC)  (Reference  5)  with  granulated  Composi¬ 
tion  B,  violent  decomposition  occurs  at  about  478  K  (205°C).  These 
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observations  led  Sewell  (Reference  6)  to  assign  the  following  values  to 
a,  A,  and  B. 

_3 

a  *  10  m/s-bar 
A  *  17.2 
B  a  0.0335/K 

These  generic  constants  give  the  infinite  burning  rate  temperature  as 
3513  X,  and  give  a  burn  rate  of  about  0.14  mm/s  for  cast  explosives  at 
3300  K.  Figure  3  illustrates  this. 


FIGURE  3.  Reciprocal  Burning  Rate  as  a  Function  of 
Explosive  Temperature. 


Pressure  Decay  Rate 


If  the  volume  under  consideration  is  vented,  the  flow  through  the 
vent  tends  to  decrease  the  pressure.  When  the  interior  pressure 
exceeds  the  outside  pressure  by  more  than  0.8  bar,  the  flow  velocity 
becomes  sonic  (Reference  7)  and  a  very  simple  expression  for  the 
pressure-decrease  results . 


dP 

dt 


a*P 


(3) 
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where 

2 

A  =  vent  area,  m 
v  * 

Cp  =  discharge  coefficient,  0.6  to  1.0 

3 

V  =  volume,  m 

a*  =  flow  velocity,  m/s 
P  =  absolute  pressure,  bars 

In  the  generic  equation,  C^,  the  discharge  coefficient  was  allowed 
to  equal  one,  i.e.,  ideal  flow.  In  actuality,  flow  through  a  square- 
edged  orifice  results  in  a  drag  coefficient  of  =0.82  because  of  the 
vena  contracta  formed  by  the  gases  exiting  the  vent  hole  (Reference  7) . 
The  sonic  flow  velocity  of  the  gases  through  the  vent  hole,  a*,  is 
computed  from  the  temperature  of  the  products,  and  is  also  affected  by 
compressible  fluid  flow.  Thus  (Reference  1) 


where,  for  a  "nominal"  combustion  gas  mixture,  with 
T  =  2500  K 
R  =  8.31434  J/mol-K 
M  =  0.028  kg/mol 
k  =  1.27, 

a*  is  approximately  725  ra/s.  Note,  that  to  do  much  better  than  the 
generic  value  requires  knowing  the  actual  composition  of  the  product's 
gas,  its  specific  heat  as  a  function  of  temperature,  and  the  actual 
flame  temperature.  These  are  different  for  each  explosive  material  and 
difficult  to  measure. 


Critical  Vent  Area 


If  the  magnitudes  of  the  pressure-decay  and  pressure-rise  terms 
are  equal,  a  critical  condition  results  in  which  the  pressure  remains 
constant.  This  condition  is  met  when  the  ratio  of  vent  area  to  burn¬ 
ing  surface  area  is  equal  to  a  constant  determined  by  the  explosive 
constants  and  the  initial  temperature. 


10 


NWC  TP  6456 


Assuming  superposition  holds,  the  pressure-rise  and  pressure-decay 
equations  can  be  combined. 

dt  =  [  (RTB  M  (A  -  BTq)  V  "  ^AvCD**^]^V  ) 

If  the  vent-area  to  burn-surface-area  ratio  is  less  than  the 
critical  value,  the  pressure  increases  exponentially;  if  greater,  the 
pressure  decreases.  Thus,  the  ratio  is  computed  as: 


SB  •  MCd.*  (A  -  BTo) 

For  the  generic  constants,  and  with  an  explosive  of  density 
1700  kg/m3,  the  predicted  critical  vent-area  to  burn-surface-area  ratio 
as  a  function  of  bulk  temperature  is  shown  for  the  Series  III  test 
condition  in  Table  3  below. 


TABLE  3.  Critical  Vent  Area  as  a  Function  of  Initial 


Explosive  Temperature  (Sg  =  11.04  in2) . 


V 

°c 

To’ 

K 

Critical  ratio, 

VSB 

Calculated  critical  vent  area 
for  VEC  test,  Av,  in2 

0 

273 

0.002161 

0.02386 

15 

288 

0.002305 

0.02545 

61 

334 

0.002896 

0.03197 

MODIFICATIONS  TO  GENERIC  EQUATIONS  FOR 
CAST  COMPOSITION  B 

The  average  burning  rate  of  Composition  B  in  over-vented  chambers 
at  ambient  temperature  appears  to  be  about  0.2  mm/s.  In  Test  VEC  1  the 
average  burning  rate  was  0.24  mm/s.  In  work  done  for  the  Defense 
Atomic  Support  Agency  (Reference  8),  the  observed  burning  rate  was 
2.5  ft/hr,  which  translates  to  0.21  mm/s.  Thus,  it  might  be  useful  to 
modify  the  [a/A  -  BT  ]  term  in  equation  (2).  Proposed  constants  for 
cast  Composition  B  are 

-3 

a  =  10  m/s-bar 
A  =  12.04 
B  =  0.0235/T 
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These  constants  give  a  burning  rate  of  0.2  mm/s  at  300  K  while  main** 
taining  the  515  K  transition  to  infinite  burning.  The  effect  of  this 
change  is  also  shown  in  Figure  3. 

A  second  change  would  be  to  use  a  more  accurate  discharge  coeffi¬ 
cient,  CL.  As  mentioned  earlier,  square-edged  orifices  at  sonic 
velocity  demonstrate  a  CQ  of  SO. 82.  This  would  be  a  good  value  to  use 
for  these  experiments,  But  for  fragment-perforated  warhead  cases,  C_ 
could  be  even  lower  (=0.6). 


Another  variable  used  in  calculating  the  critical  vent  area  is  the 
burning  surface  area,  SR.  For  a  system  burning  linearly  (cigarette 
fashion),  Sg  is  constant  if  ignited  over  the  whole  surface  instanta¬ 
neously.  However,  for  the  ignition  system  used  in  these  experiments, 
only  a  strip  of  the  cylindrical  explosive  surface  was  initially 
ignited.  Here,  Sg  should  be  calculated  as  a  function  of  time.  In  the 
worst  case,  the  explosive  could  burn  at  its  linear  burning  rate  in  the 
X,  Y,  and  Z  directions  simultaneously.  For  this  particular  geometry, 
Sg  as  a  function  of  time  was  calculated  for  the  over-vented  situation. 
Initially,  the  burning  surface  is  smaller  than  the  total  exposed  explo¬ 
sive  surface  of  11.04  in2.  At  96  seconds,  the  burning  surface  is 
equivalent  to  11.04  in2  (7126  mm2,  the  total  exposed  surface).  At 
155  seconds,  the  trough-shaped  burning  surface  reaches  a  maximum  value 
of  8705  mm2  and  by  280  seconds  has  returned  to  7126  mm2.  Figures  4 
and  5  illustrate  this.  Note  that  if  central  point  ignition  had  been 
used,  there  would  have  been  a  170-second  delay  to  S0  =  7126  mm2  and  an 
overshoot  of  almost  87%  at  230  seconds.  (The  burning  surface  would  be 
spherically  shaped.) 


The  effects  of  changing  the  discharge  coefficient  in  the  generic 
equation  from  1.0  to  0.82  is  shown  in  Figure  6.  Note  that  VEC  data 
points  for  Test  No.  2  and  5,  which  rapidly  pressurized,  fall  just 
beneath  this  line.  Test  No.  3  is  considered  anomalous,  as  stated  in 
Table  2.  Applying  the  proposed  modifications  to  the  burning  rate  equa¬ 
tion  for  Composition  B  shifts  the  boundary  between  violent  pressure 
rupture  and  quiescent  burning  upward  to  larger  values  of  A  /SR.  This 
appears  to  be  a  slightly  conservative  representation  of  tKe  boundary 
from  a  safety  point  of  view  since  all  rapid  burns  lie  below  it  and  the 
quiescent  burns  lie  upon  it.  This  is  in  qualitative  agreement  with  the 
Air  Force  work  (Reference  8)  shown  in  Figure  7.  They  found  quiescent 
burns  in  the  range  of  0.121  to  0.003  A  /SR,  while  violent  pressuriza¬ 
tion  occurred  in  many  tests  in  the  range  of° 0.0058  to  0.0028  Av/SR. 


THERMAL  MODELING  OF  HEAT  FLOW 

As  was  seen  in  the  previous  sections,  the  required  critical 
vent  area  to  burning  surface  area  increases  as  the  bulk  explosive 
temperature  increases.  One  can  visualize  an  impacted  explosive  round 
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FIGURE  4.  Burning  Surface  Area  of  Composition  B  Versus 
Time  for  Various  Ignition  Source  Geometries. 
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FIGURE  5.  Shape  of  Burning  Surface  as  a 
Function  of  Time  (X-Z  Plane). 


*******^ 


NWC  TP  6456 


0.006 

0.004 

£ 

^  0.003 
0.002 

0.001 


0L 


c* 

z 


0.06 


0.06 


0.04 


0.03 


h  > 


0.02  h 


0.01  h 


CRITICAL  CONDITION  - 

COMPOSITION  B  CONSTANTS  -  Cn  -  0.82 

/ 


CRITICAL  CONDITION  - 
3ENERIC  CONSTANTS  -  C0  -  0.82 
(NO.  6,  4.4  •) 

CRITICAL  CONDITION  - 
GENERIC  CONSTANTS  -  C0  -  1 
(NO.  2,  7.6  ft) _ 

O  VEC  TESTS  ~ 

QUIESCENT  BURNS 
•  VEC  TESTS  - 

RAPID  PRESSURIZATION 

(NUM8ERS  IN  PARENTHESIS  ARE 
VEC  TEST  NO.  AND  BURN  TIME) 


250  300  350  400 

INITIAL  EXPLOSIVE  TEMPERATURE  (T0).  K 


FIGURE  6.  VEC  Test  Results  Compared  to  Predicted  Critical 
Vent  Area  for  Sg  =  11.04  in2. 


that  ignites  and  burns  in  the  over-vented  condition  for  some  time;  but 
because  of  heat  transfer  down  the  wall  of  the  device  and  radiation 
heating  of  the  bulk  explosive,  the  explosive  temperature  rises,  the 
reaction  rate  increases,  and  the  available  vent  is  suddenly  too  small 
to  prevent  pressure  rise. 


To  determine  the  magnitude  of  this  effect,  Stroebel  and  Graham 
(Reference  9)  (see  Appendix  B)  performed  a  simplified  thermal  analysis 
of  the  VEC  tests.  In  this  initial  work,  the  linear  burning  rate  of 
Composition  B  was  assumed  to  be  0.1  mm/s.  Thermal  contours  were  gener¬ 
ated  using  the  SINDA  computer  program.  These  contours  show  definite 
preheating  of  the  bulk  explosive  ahead  of  the  flame  front.  While  the 
burning  rate  was  too  slow,  and  a  kinetic  equation  for  rate  of  change 
of  burning  velocity  with  temperature  was  not  incorporated  into  the 
routine,  it  is  significant  that  the  steel  wall  was  heated  ahead  of  the 
flame  front  and  that  the  explosive  near  this  wall  was  also  signifi¬ 
cantly  preheated.  Figures  8,  9,  and  10  show  the  results  of  these  cal¬ 
culations  at  three  different  response  times. 
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CONCLUSIONS 


In  response  to  the  questions  posed  in  the  introduction  of  this 
report: 

1.  Two  relationships  between  critical  vent  area  and  burning  sur¬ 
face  area  have  been  developed — one  for  generic,  cast  explosives  and  one 
specifically  for  Composition  B. 

2.  The  critical  vent-area  to  burning-surface-area  ratio  (A  /SR) 
is  affected  by  the  bulk  temperature  of  the  explosive.  The  higher  the 
bulk  temperature,  the  faster  the  burning  rate  and  the  larger  the  vent 
required  to  prevent  rapid  pressurization,  as  indicated  in  Equation  (2). 

3.  Delayed  violent  reactions  have  been  observed  in  burning, 
veated  munitions  and  in  small-scale  tests.  It  appears  that  the  change 
in  critical  vent  size  as  the  burning  surface  increases  in  area  and  in 
temperature  plays  a  key  role  in  these  delayed  reactions. 
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FIGURE  8.  Cross  Section  of  Steel  Cylinder 
Showing  Isotherms  in  Burning  Explosive  at 
60  Seconds. 
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UJ 


FIGURE  9.  Cross  Section  of  Steel  Cylinder 
Showing  Isotherms  in  Burning  Explosive  at 
180  Seconds. 
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RADIAL  POSITION,  INCHES 


FIGURE  10.  Cross  Section  of  Steel 
Cylinder  Showing  Isotherms  in  Burn¬ 
ing  Explosive  at  300  Seconds. 


4.  It  is  apparent  that  either  the  generic  venting  equation 
with  CD  =  0.82  or  the  modified  equation  (also  with  C^  =  0.82)  would 
adequately  represent  the  boundary  between  quiescent,  vented  burning; 
and  violent,  undervented  pressure  rupture.  In  light  of  the  Air  Force 
work  cited,  and  due  consideration  for  a  margin  of  safety,  an  Av/Sg 
ratio  of  >0.006  at  an  explosive  temperature  of  300  K  should  prevent 
violent  pressurization  in  Composition  B-loaded  rounds. 
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Appendix  A 

SAMPLE  TEST  DATA  FOR  VEC 
SERIES  III  TESTS 


FIGURE  A-l.  VEC  Test  1,  Pressure-Time  History. 

A  *  0.0430  in2,  T  =  15°C. 
v 


♦ 


FIGURE  A-2.  VEC  Test  2,  Pressure-Time  History. 

A  =  0.0287  in2,  T  =  15°C. 
v 
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THERMAL  MODELING  OF  VEC  TESTS  USING  THE  SINDA  PROGRAM 
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DEPARTMENT  OF  THE  NAVY 
NAVAL  WEAPONS  CENTER 
CHINA  LAKE,  CAUH'CIINIA  »»»*J 


in  RtPLV  mrcn  »a. 

3242/FAS: roc 
Keg  3242-27-83 
3  Mar  1983 


From:  F.  A.  Sttobel,  Thermal/Structurus  Branch  (Code  3242) 

To:  K.  J.  Graham,  Warhead  Dynamics  Branch  (Code  3835) 

Via:  Head,  Thcrmal/Structures  Branch  (Code  3242)  AW/ 

Subj:  Composition  B  thermal  analysis 

Ref:  (a)  TRW  report  11027-6003-R0-0O,  dtd  Sept  1980. 

(b)  Chapman,  "Heat  Transfer,"  MacMillan  Publishing  Co.,  1974. 

(c)  Aerotherm  report  UM-75-68,  dtd  Dec  1975. 

Enel:  (1)  Sketch  of  test  setup 

(2)  Sketch  of  conduction  nodal  network 

(3)  Material  properties 

(4)  Explosive  thermal  response,  time  ■  300  seconds 

1.  The  purpose  of  this  memorandum  is  to  describe  the  methodology  used  in  con¬ 
ducting  the  thermal  analysis  for  the  Composition  B  explosive  tests  conducted 
by  Code  3835.  The  analysis  was  used  to  calculate  temperature  response  of  the 
Composition  B  explosive  due  to  heat  conduction  from  the  hot  surface  and  through 
the  steel  walls  of  the  holder.  This  analysis  was  completed  in  November  1982 
and  the  results  were  forwarded  to  Code  3835. 


2.  The  temperature  calculations  were  made  with  the  use  of  the  Systems  Improved 
Numerical  Differencing  Analyzer  (SINDA,  reference  (a)).  Enclosure  (1)  shows  a 
sketch  of  the  test  setup  and  the  phenomena  considered  in  the  analysis.  Besides 
heat  conduction  throughout  the  explosive,  liner,  and  steel,  the  thermal  model 
included  convective  heating  of  the  inner  surface  of  the  steel  cylinder  exposed 
to  the  hot  combustion  gases.  In  addition,  the  model  included  cooling  of  the 
steel  cylinder  due  to  radiation  and  natural  convection  to  the  surroundings. 

3.  Since  the  part  being  modeled  is  axi-symmetric ,  a  two  dimensional  conduction 
model  was  sufficient.  Because  the  explosive  was  burning  at  its  upper  surface 
and  receding  at  the  rate  of  0.1  mm/r.ec,  it  was  necessary  to  continuously  regrid 
the  conduction  model.  A  sketch  of  the  nodal  network  for  the  conduction  model 
is  shown  in  enclosure  (2).  The  Comp  B,  J.incr,  and  steel  were  comprised  of  72, 
28,  and  38  nodes,  respectively.  An  extra  fine  grid  was  applied  to  the  top  sur¬ 
face  and  outer  edges  of  the  explosive  because  these  were  the  areas  where  larg¬ 
est  temperature  gradients  would  occur.  In  order  to  keep  the  nodes  aligned  all 
three  parts  had  to  be  continuously  regrldod.  Since  the  node  centers  changed 
after  each  time  step,  an  interpolation  scheme  had  to  be  used  to  find  the  temper¬ 
atures  at  the  new  node  center  locations. 


4.  Convective  heating  rates  to  the  inner  wall  of  the  steel  cylinder  were  cal¬ 
culated  from  a  method  for  pipe  flow  described  in  reference  (b).  The  accuracy 
of  this  method  is  questionable  and  this  could  be  an  area  of  improvement  if 
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3242/FAS:cec 
Reg  3242-27-83 


Subj:  Composition  B  thermal  analysis 

*•  additional  analyses  were  needed.  To  calculate  the  convective  heating  rates 

it  was  necessary  to  know  the  properties  of  the  hot  combustion  gases.  These 
properties  were  obtained  through  use  of  the  Aerotherm  Chemical  Equilibrium 
#  computer  code  (ACE,  reference  (c)). 

5.  The  thermal  property  data  used  for  the  analysis  is  contained  in  enclo¬ 
sure  (3).  More  detailed  property  data  for  the  Comp  B  and  liner  materials 
would  add  to  the  accuracy  of  the  calculations.  In  addition,  modelling 
the  decomposition  and  burning  of  the  liner  material  would  also  add  to  the 
accuracy  of  the  calculations. 

6.  In  conclusion,  a  reasonably  good  thermal  model  for  the  Comp  B  explosive 
tests  has  been  developed.  This  model  takes  into  account  the  most  important 
phenomena  occurring.  A  typical  plot  of  results  is  shown  in  enclosure  (4). 
This  plot  shows  steel  temperatures  and  isotherms  in  the  explosive  after  a 
burn  time  of  300  seconds.  Possible  areas  for  improvement  in  the  model  have 
been  outlined  in  the  preceding  paragraphs. 

^■d.  Jfiubd. 

P.  A.  Strobel 


Copy  to! 
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MATERIAL  PROPERTIES 


4 


Steel  Temperature  (#F) 


31.4 

166.4 

391.4 

751.4 
1111.0 
1291.0 

Temperature  (8F) 


32.0 

392.0 

572.0 

752.0 

1112.0 

1472.0 

1832.0 

2192.0 

Composition  B  Temperature  (°F) 


30.0 

86.0 

122.0 

158.0 

212.0 

5000.0 


Density  X  Specific  heat 
(Btu/in3-*F) 

0.2982 

0.03408 

0.03834 

0.0426 

0.0482, 

0.05679 

Conductivity  (Btu/in-sec-°F  x  10**) 


4.861 

4.861 

4.630 

4.398 

4.167 

3.704 

3.704 

3.935 

Density  X  Specific  heat 
(Btu/in^-’F) 

0.0186 

0.0186 

0.0191 

0.0203 

0.0207 

0.0207 


4 


Conductivity  for  all  temperatures  ««  3.484  X  i0-&Btu/in-sec- WF) 


Liner 


Density  X  Specific  heat  ■  0.0158  Btu/in3-°F 
Conductivity  ■  9.9  X  10*”6Btu/in-sec-°F 


V 

J 


All  temperatures 


Enclosure  (3) 
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